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Quantum master equations are an important tool in quantum optics and quantum information 
theory. For systems comprising a small to medium number of atoms (or qubits), the non-truncated 
equations are usually solved numerically. In this paper, we present a group-theoretical superoper- 
ator method that helps solving these equations. To do so, we exploit the S(7(4)-symmetry of the 
respective Lindblad operator and construct basis states that generalize the well-known Dicke states. 
This allows us to solve various problems analytically and to considerably reduce the complexity of 
problems that can only be solved numerically. Finally, we present three examples that illustrate the 
proposed method. 
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I. INTRODUCTION 



Quantum master equations have been used to address a variety of problems from quantum optics, such as the 
description of laser systems and the phenomenon of resonance fluorescence [H EJ 0] . In quantum information theory, 
they are used to study the decoherence of entangled quantum states [7J. For a finite number Z of 2-level atoms (or 
qubits), quantum master equations involve non- unitary terms of the Lindblad form: 
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Here P is the reduced density operator of the system under consideration, and a± and a^ are the usual Pauli 
matrices, acting on atom i whose states are |1 > or |0 >, representing, respectively, the atom being in the excited 
state or in the ground state. B and C are decay constants and s is the pumping parameter pj]. It varies from s = 
for pure damping to s = 1 for full laser action. 

To solve the corresponding quantum master equations, three approaches have been taken: First, one focuses on 
the case of one atom. Second, one truncates eq. ([I]) and derives semi-classical models. Third, one employs numerical 
simulation methods such as the quantum trajectory method [4|. While the latter method is very popular, it should 
be noted that the numerical complexity of the problem increases exponentially with the number of atoms, and so 
numerical methods soon become unfeasible. To proceed also in these cases, we propose an analytical approach that 
exploits the symmetry of the Lindblad operator. 

We will show in this paper that Lindblad operators such as the one in eq. ([1]) are invariant under SU (4) transforma- 
tions. This observation allows us to construct appropriate basis states for the density operator P. These basis states 
can be classified according to their symmetry type Y, which is preserved by the unitary and non-unitary dynamics 
of the system. Our approach has two advantages: First, being able to express the Lindblad operator in terms of 
the generators of SU(4) allows us to solve a large class of problems analytically. Second, the observation that the 
dynamics preserves the symmetry type Y leads to a considerable reduction of the number of required basis states for 
num eric al calculations. For example, for fully symmetrical states, the number of basis states is of the order Z 3 (see 



eq. (39)), which has to be compared to 4 if the symmetry of the dynamics is not taken into account. 
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The remainder of this paper is organized as follows. Sec. [TT] introduces Dicke states and shows that they do not 
suffice to solve quantum master equations that involve non-unitary terms such as those in eq. . Sec. Ill introduces 
generalized Dicke states and explicitly constructs fully symmetrical basis states. Sec. |IV| shows how the proposed 
method can be used to solve quantum master equations. Sec. [V] present three examples that illustrate the proposed 
method. Finally, Sec. 



VI suggests a number of further applications that we plan to address in future work. 



II. DICKE STATES 



As will become apparent below, the new basis states that we propose in this paper generalize the well-known Dicke 
states [5]. To show the need for a generalization, we first introduce Dicke states and show that they do not suffice 
to solve quantum master equations involving terms such as those occurring in eq. 0. We consider a system of Z 
2- level atoms, of which n + atoms are in state |1 > and n_ = Z — n + atoms are in state |0 >. Using the analogy to 
the standard SU(2) angular momentum eigenstates, Dicke states are characterized by two quantum numbers, S and 
S3. They are defined as follows: 

\S,S 3 >=S\{n+,n_}> (2) 

Here S is the symmetrization operator, S = Z/2 and n± = S ± S3 with —S < S3 < S . As an example for a Dicke 
state of three atoms, we calculate 

|3/2, 1/2 > = S|{2,1}> 

= — (|110 > +|101 > +|011 >). 
v3 

Dicke states are eigenstates of S 2 (with eigenvalue S(S + 1)) and S3 (with eigenvalue S3). Raising and lowering 
operators (S±) and S3 can be defined accordingly: 
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They satisfy the commutation relations 

[S+,S_]=2S 3 , [S 3 ,S±]=±S±. (4) 

Using Dicke states, it is straight forward to construct basis states for density operators for a system composed of 
Z atoms: 

P^fJ = \Z/2, M >< Z/2, M'\, for M, M' = -Z/2, ...,Z/2 (5) 

Once the number of atoms is fixed, these states are characterized by two quantum numbers, M and M'. The basis 
has the dimension 

D D {Z) = {Z + lf. (6) 

Unfortunately Dicke states do not suffice to diagonalize quantum master equations for Z atoms that involve non- 
unitary terms such as those in eq. ([l}. To show this, we first note that the first two terms in eq. ([T]) work 

fine as J2i=i °+ = Z/2 + S3 and J2i=i 0- 0+ = Z/2 — S3. However, calculating Y^=i 0_ -P0+ for, say, 
P 3 ( aJ/2 = I 3 / 2 ' 3 / 2 >< 3/2, 3/2j = |111 >< 111| yields 

3 

£ a W P^ 3 3 3 ) 2 = |011 >< Oil] + |101 >< 101] + |110 >< 110|, (7) 

t=i 

which cannot be expressed as a superposition of basis states Pp^ M' ■ 

If one nevertheless wants to use Dicke states to tackle such problems, one can truncate eq. (fTl) and use the operators 
S± and S3 to arrive at an expression L of the Lindblad form that respects the condition TryLP) = 0. For C = B/2, 
i.e. neglecting the dephasing term, one obtains: 

L a iZ) P = ---(l-s) {S+S-P + PS+S- - 2 S-PS+) - -s (S_S+P + PS-S+ - 2 S+PS-) (8) 
_ z 
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(Note that the natural replacement of the dephasing term by P — S3PS3 is not possible as S3 7^ 1 for Z > 2.) We 



conclude that whether or not eq. ([8 1 is a good approximation to eq. ([!]) (for C — B/2) has to be decided on 



case-by-case basis by comparison wit 
prove useful 



a 



1 the exact solution. To find the exact solution, generalized Dicke states will 



III. GENERALIZED DICKE STATES 

To begin with, let us study the algebraic properties of superoperators such as those occuring in eq. 0. We will 
see that these superoperators satisfy the commutation relations of the generators of the special unitary group SU(A). 
This observation will help us to construct appropriate basis states. These basis states are classified according to their 
symmetry type Y , which is typically represented by a Young tableau. 

A. SU(4:) Superoperators 

We define the following superoperators: 

z „ z 
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M ± P:=j:^p l -±^l , M,P:=\j:^P l -±^- (9) 
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15 of these 18 superoperators are linearly independent. In fact, it is easy to see that A/3 = Q3 + S3 — M.3, U3 = 
—S3 + M3 and V3 = Q3 — .M3. The 15 independent superoperators are the generators of the Lie group 577(4) 
(see [SI [§])■ Table II in Appendix A lists all commutation relations between the superoperators. Here we focus 
on some important commutation relations that will help us to construct the basis states. To simplify notation, let 
O := {Q,Y,,M,AT,U,V}. We then note that for all X G O: 

[X+,X-] = 2X3 (10) 
[X 3 ,X±] = ±X ± (11) 

That is, we have identified six SU (2) subgroups. Note also that 

[Qi,Ej] = [MuAfj] = [U t ,V 3 ] = G {±,3}, (12) 

i.e. the sw(2)-subalgebras for Q and X etc. are "orthogonal". Next, it is useful to study the quadratic superoperators 

X 2 = X-X + +Xl+X 3 , (13) 

for all X G O. These superoperators satisfy the commutation relations 

[A 2 ,A' 2 ]=0 and [A 2 ,A 3 ]=0, (14) 

for all X,X' G O. Moreover, we note that 

[X 3 ,X' 3 }=0, (15) 
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TABLE I: The action of the SU(4) superoperators on the elements of the fundamental representation. 



for all X, X' e O. Hence, there are joint eigenfunctions of Q 2 , Q 3 , E 2 , S3, M 2 , M 3 , TV 2 , A/3, U 2 , U 3 , V 2 and V 3 . 
Before constructing these eigenfunctions (or basis states) , we note some important duality relations that hold between 
our superoperators. The dual conjugate L to a superoperator L is defined as 

Tr{0(LP)} = Tr{{LO)P] (16) 

for all states P and all observables O (5]. It then follows from the definition of our superoperators that 

Q± = Q=f , 2s = Qs 

s± = S T , S 3 = —S 3 

M ± =U T , M 3 =U 3 (17) 

-A4 = V T , Af 3 = V 3 . 



B. Basis States 



The fundamental representation of the group 5/7(4), adapted to the present case, is explicitly given by 

u := |1 >< 1| , d := |0 >< 0| 

s := |1 >< 0| , c := |0 >< 1|. (18) 

The names allude to the SU(A) quark model. (As Georgi notes, the Flavor SU(4) model was actually never very 
successful, but this is not relevant for the present discussion, see [B].) Table I shows how our superoperators (for Z = 1) 
act on these states. All higher order representations can then be obtained from the fundamental representation and 
the symmetry type Y, which is conveniently characterized by a Young tableau comprising four rows of pi,p2,p 3 and 
P4 boxes with pi + p 2 + p 3 + P4 = Z. 



It is important to note that the symmetry type Y of a state is preserved under the action of the superoperators 
defined in (|9j): if one of the superoperators (or a linear combination thereof) is applied to a state of symmetry type 
Y, then the resulting state also has the symmetry type Y. Hence, the non-unitary dynamics preserves the symmetry 
type Y as the terms on the RHS of eq. |T]) can be expressed in terms of the superoperators defined in (J9J) . The 
same holds for the unitary dynamics. To see this we note that the atomic part of the Hamiltonian H involves terms 
proportional to S± and S 3 . The commutator [H,P], which captures the unitary dynamics in the quantum master 
equation, then involves terms proportional to S± P = (M± + Af±) P, P S± = (U± + V±)P, S 3 P = (Q 3 + S 3 ) P and 

PS 3 = (Q 3 -S 3 )P. 

Let us now identify the basis states. For the configuration u a d 13 s 1 c s (with a + /3 + 'f + 5 = Z), these are given by 

P Y =S Y (u a d^s->c s ), (19) 

where the symmetrizer S Y makes sure that the state has symmetry type Y. S Y is the product of symmetrizers and 
anti-symmetrizers which can be read-off from the corresponding Young tableau. (Note that the symmetrizers have to 
be applied first, then the antisymmetrizers, so that the resulting state is antisymmetrical in the corresponding indices, 
but not necessarily symmetrical in all indices which were symmetrized [9].) 

To construct the basis states, one can also start with a state of highest weight (for a given symmetry type Y) and 
apply the raising and lowering operators Q± etc. to generate all other states. For full symmetry (i.e. for the Young 
tableau (Z, 0, 0, 0)), the states of highest weight are u z ,d z , s z and c z . 
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follows from two considerations. First, the explicit form of the fundamental representation (see eq. (18)) implies that 

Tr(u) = Tr{d) = 1 and Tr(s) = Tr(c) = 0. (20) 



Hence, only composite states containing it's and d's have a non-vanishing trace. Second, because of the previous 
consideration, Young tableaux for non-fully-symmetrical states have exactly two rows. (A third or forth row would 
require elementary s or c states to fulfill the symmetry requirement.) Hence, each term in the corresponding state 
contains at least one antisymmetrized u/d pair, and so the trace of such states vanishes. 

In the remainder of this paper, we focus on fully symmetrical states and leave a detailed analysis of states of mixed 
symmetry and fully antisymmetrical states to a sequel to this paper. We reserve the name generalized Dicke states 
for fully symmetrical basis states of trace 1. (For a different way to generalize Dicke states, see [8].) Strictly speaking, 
we should not use the word "state" at all if the trace vanishes. For convenience, however, we will do so anyway. This 
also makes sense as states of vanishing trace can be created by the dynamics as admixtures to states of trace 1. 



by 



C. Fully Symmetrical States 

For ease of notation, we denote the symmetrizer for fully symmetrical states by 5, so that the basis states are given 



p( s) = S(u a dP s J c s ). 



(21) 



To specify suitable quantum numbers that label the basis states, we note that S commutes with all our superoperators. 
Using Table I, we then derive: 
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With eq. (13) and requesting (or defining), as usual, that X 2 



s) 



I £ O, we obtain from the first line of eqs. (22) 



(x + 1) P^ and X 3 P( 
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Hence (neglecting the negative solution) 
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Similarly, we obtain: 
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From this we conclude that each term in a fully symmetrical basis state has (i) the same numbers of l's (= 2m) and 
O's (= 2n) on the RHS (i.e. the < -|'s) of the density operator, (ii) the same numbers of l's (= 2u) and O's (= 2v) on 
the LHS (i.e. the |- >'s) of the density operator, and (hi) the same number of — or 1 — 1 coincidences (= 2q) and 
1 — or — 1 non-coincidences (= 2a) on the LHS and the RHS of each term in a basis state. These observations 
give our quantum numbers some intuitive meaning. 

All fully symmetrical basis states for a given number Z of atoms can then be characterized by three quantum 
numbers. We choose q, q 3 and 03 from which we can calculate a, /3, 7 and 6 in turn: 

a = q + q 3 , /? = <?- qz 

7 = a + a 3 , 8 = a - cr 3 , (28) 

with a = Z/2 — q. Let us now explore the range of the values of these quantum numbers. From [Q3, Q±] = ±Q± we 
conclude that Q+ raises q 3 by one unit until ft = 0, i.e. until q 3 = q. Similarly Q_ lowers q 3 by one unit until a = 0, 
i.e. until q 3 = —q. Hence, q 3 ranges from —q to q. Similarly, cr 3 ranges from —a to a. From eq. (24) we conclude 



( z) 

• : q ranges from to Z/2 in steps of 1/2 unit. We denote the basis states by Pq,qi,a 3 and arrive at the following 
list of the three quantum numbers that characterize them: 



1 

q 3 -Q 3 P. 
S3 P, 
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9,93,0-3 
(Z) 



DP 
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y3 -"9,93, 0-3 ! 
?(Z) 



?,93,<X3) with 9 = 0, 1/2, . 
with —q < q 3 < q 



.,Z/2 



q,q 3 ,a 3 — C 3 Pq^qi,a 3 , with —a < O3 < (J and (7 = Z/2 - q 

Note that the quantum numbers m, 7713 etc. can be expressed in terms of q, q 3 and 03. Given the symmetry of our 
approach, we could have also used the eigenvalues of the M./M or U /V superoperators to label our basis states. The 
Q/S labeling is, however, most convenient for the applications we have in mind. 

Generalized Dicke states are fully symmetrical basis states of trace 1. For them, 7 = 5 = and hence a + j3 = Z. 
Using eq. (24), we conclude that generalized Dicke states are characterized by the quantum numbers q = Z/2, q 3 — 
—q, . . . , q and 03 = 0. All other fully symmetrical basis states have a vanishing trace. To see how the basis states 
look like, we calculate a generalized Dicke state with highest weight, 
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Applying Q_ to Pi^ 0/on , we obtain |011 >< 011 



su 2 ) 

o 

111 > (< iioj+ < 101 l+ < 011 D- 

110 >< 110 1 (see eq. M), which equals 



3 P. 
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3/2,3/2,0' vvc vu^aiii |<jj.i UJ-J-I I" 1 101 >< 101 1 

3/2 1/2 0" us now study the action of the raising and lowering operators on our basis states more generally. For 
Q± and £±, we obtain: 



Q± P {z) 



^±P ( q 
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93,03 



= (?T <73) Pq Z q\±l^ 
( \ n(Z) 



(30) 
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To calculate the action of A4 + on our basis states, we observe that 

M+ S (u a <tf s 7 c 5 ) = 5S(u a+1 d p s 7 c 5 - 1 ) =: 6S(u a ' di 3 ' s 7 ' c 5 '). 



(31) 



Hence, with eqs. (28), q' = q + 1/2, q' 3 = q 3 + 1/2 and a' 3 = (T 3 + 1/2. Similarly for the other superoperators, so that 
we finally obtain: 



M ± P ^ 3 = ( m ^ m ^ P g±i/2, q3 ±l/2,, 3 ±l/2 
N±Pq Z q 3 ,cr 3 = ( n T "3) Pq^l/2,q 3 ±l/2,a 3 ±l/2 



(32) 



9Tl/2,93±l/2,o 3 =Fl/2 

= (vTv 3 )P { q Z ± \ /2iq3±1/M/2 



The expressions for the new values of q$ and 03 after a certain raising or lowering operator acted on the corresponding 
state can also be obtained from commutation relations such as [S3,A^±] = ±1/2A / J± and [£3,.M±] = ±1/2 A4±. 
However, calculating in this wa y how the value of q changes is more difficult as Q 2 is a quadratic operator. But a 
trick helps. Inspired by eq. (24), we define the superoperator Q by 



G:= \ (Z+Q33) 



with 



(33) 



Q33 P := 4 ] P4 } = (4^3 - 2 (Q 3 + S 3 )) P. 
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It is easy to see that 



and that, for fully symmetrical states, 



[Q, Qi] =0 for i G {±, 3} 



Si 3 , 



(2) _ 



9,93,0-3 



y 9,93,03 
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Besides, [Q, A4±] = ±1/2 A4±. Hence, Ai± raises/lowers the value of the quantum number q by 1/2 unit for fully 
symmetrical states. Similarly, linearized superoperators X (for all IgO) can be constructed in a straight forward 
way from eqs. (26). Using these superoperators simplifies calculations considerably, but note that they only work for 
fully symmetrical states. 

Next, we define the dual states, 



p(Z) = p{Z) 

9,93,03 1 9,93,-03' 



which satisfy the biorthogonality relation 



Z) p(Z) 



S q q> 5 q3 q ^ 5 a3 CT £ . 



(37) 



(38) 



To conclude this section, we determine the dimension of the space containing all fully symmetrical states. By simple 
counting or using the hook rule for Young tableaux, we obtain 



D GD (Z) = i(Z+l)(Z + 2)(Z + 3). 
6 



(39) 



That is, Dgd(5) = 56,Dgd(10) = 286 and Dqd(2Q) = 1771. These numbers have to be compared to A z if one 
proceeds by brute force. Note, for example, that 4 10 w 3670 x 286, i.e. for ten atoms, we save already a factor of 
3670. Compared to the number of Dicke states for the same number of atoms (see eq. (|6)), we obtain, for example, 
Dgd(10)/-D,d(10) « 2.4 and Dgd(50) / Dd(50) « 9, which indicates only a relatively moderate increase. Hence, our 
basis states will be useful for numerical solutions of quantum master equations. In the next section we show that 
some quantum master equations can even be solved analytically. 
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IV. SOLVING QUANTUM MASTER EQUATIONS 

Consider first a quantum master equation with a Lindblad operator given by eq. ([TJ for C = B/2 (i.e. without the 
dephasing term). That is, we want to solve the equation 



dP 



L (Z)p 



B 



-(1 - s) KV W P + Pafa { y - 2aWPaf) 



i=l 



B 



i=l 



Using the superoperators defined in ([9]), this equation can be compactly written as 

1 dP 

-— = [-Z/2 + (1 - s) Q- - (1 - 2s) Q 3 + s Q + ] P. 



We introduce the dimensionless time r := Bt and solve eq. (41 1 formally: 



P(r) 



_ p -Z/2t Ml-s) Q_-(l-2s) Q 3 +s Q+)t 



P(0) 



(40) 



(41) 



(42) 



Here P(0) is the initial state of the system. We now use the Baker-Campbell-Hausdorff formula to factorize the second 
exponential and obtain the following two formulas (which one is most convenient will depend on the initial state). 



P(r) = e -2/^^(r)Q +e B s WQ 3e C l( r)g_ F(fl) 



(43) 
(44) 



with 



A s {t 
B s {t 
C s (t 

D s {t 
E s {t 
F s (r 



sf(r) 



l-s/(r) 

-r-21og(l- S /(r)) 
(l-s)/(r) 
1-s/(t) 
(l-s)f(r) 



= Ai_,(t) 

0/W) = 

Ci_.(t) 



and 



l-(l- s )/(r) 

r + 2 log(l - (1 - s) /(r)) = -Px_ s (t) 
l-(l- S )/(r) 



f(r) 



(45) 



(46) 



Master equations involving other or more superoperators from (|9| can be formally solved in a similar way. Here we 
only show how the dephasing term from eq. ([T]) can be included. To do so, we first note that 

i=l 

= 5(1 - 2C) (Z/2 -Q)P 

=: B{l-2C)tP, (47) 



with C :— C/B. Next, we recall that Q (and hence E) commutes with Q± and Q3 (see eq. (35)). Hence, the solution 
of the quantum master equation including the dephasing term simply adds the factor 



P dep h(r) = e^( z / 2 -^ 



(48) 
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to the right hand sides of eqs. (43 1 and (44). Note that generalized Dicke states (for which q — Z/2) are not affected by 
the dephasing term. In the following illustration and the applications in the next section, we disregard the dephasing 
term. 

To illustrate the method developed in this section, let us calculate the time evolution of the fully symmetrical state 



± q,q 3 ,cr 3 under the influence of pure decoherence, i.e. for s — 0. Using eq. (43), the resulting state is given by 

p (T)=e -Z/2r e -rQ 3 e /Wfl-^ iffji 

and after some algebra, we obtain the compact expression 



(49) 



P(r) 




/ fe (r)(l-/(r)p- fc P ( 



(Z) 

g,q 3 —k,a 3 



(50) 



Using eq. (50), we can easily calculate the time evolution of the atomic inversion < 5,3 >=< Q3 >= Tr(Q 3 P(t)) of 



(Z) 

a decaying fully symmetrical state. For the initial state ( w ^ n even we obtain < S3 > jZ = 1/2 (e~ T — 1). 

1/2. It is instructive to compare the 



For the initial state P, 



Z/2,Z/2fl 



Z/2 

= |1 . . . 1 >< 1 . . . 1|, we obtain < S 3 > jZ = 



latter result with the expression that one obtains if one uses Dicke states and replaces the present Lindblad operator 



by the one given in eq. 



for s = 0. In the Dicke basis, the initial state is given by P, 



inversion for Z = 2 is < S3 > /2 = (1 +r) e 



-2t 



(D,Z) 
Z/2, Z/2 



and the resulting atomic 



1/2, which suggests a much faster decay. 



V. APPLICATIONS 



Let us now consider three applications of our formalism. We first study the propagation of a Bell state for an 
arbitrary pumping parameter s (Sec. VA). Then we calculate the decay of a GHZ state (Sec. VB). Finally, we 
calculate the eigenvalues of the Lindblad operator (Sec. VC). 



A. The Propagation of a Bell State 



We calculate the propagation of the Bell state \ 1 i> + >: 



1 01 > +|10 >) for arbitrary s. To do so, we first note 



that P B (0) := |*+ >< *+| = S(ud - 



p|q o + fp'o an d obtain from eq. (50) for the time evolution of Pb(0): 



j(2) 



Pb(t) = 61 (r) p£ + b 2 (r) P« + b 3 (r) P™ + b 4 (r) P& 



(51) 



The weights are given by h(r) = s f(r) [1 - (1 - s) /(r)], 6 2 (r) = 1 - /(r) [1 - 2s(l - s) /(r)], b 3 (r) = (1 - s) /(r) [1 - 
s/(r)] and 6 4 (r) = 1 - /(r). For s = 1/2, P b (t) simplifies to 



Pb{t) = \(l- er^) P[% + \ (1 + e--) Pfl, + 1(1- e-*) P[% fi + e- < 0) 



(52) 



which converges asymptotically to 1/4 P^l + 1/2 P^ + 1/4 Pfl_ 1 . More generally, Pb (t) converges asymptotically 
to 



,(2) 



,(2) 



P B (r) = S 2 Pg + 2 S (1 - 8 ) Pgo + (1 - .s) 2 Pl% t , 



(53) 



Let us now calculate the von Neumann entropy of Pb(t) as a function of time r for different values of the pumping 
parameter s. To do so, we diagonalize the density matrix from eq. (51 1 and calculate its eigenvalues for i = 1, . . . ,4. 



The von Neumann entropy is then given by S_b(t) = —Tr{PB{r) log Pb(t)) = — J2i=i ^ l°g^t- (The base of the log 
is 2.) It is depicted in Figure 1. 



B. The Decay of a GHZ State 



Next, we consider the case of pure damping (i.e. s — 0), and propagate a GHZ state through time. To do so, we 
express the density operator corresponding to the GHZ state 



|GHZ >= -p (| 111 > -|000 >) 
v2 



(54) 
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2.5 

2 



1.5 




1 2 3 4 5 



FIG. 1: The von Neumann entropy Sb(t) of Pb(j) as a function of time r for s = .5 (dotted), s = .1 (dashed) and s = .9 
(solid). 



in our fully symmetrical basis for Z — 3. We obtain: 
Pghz 



|GHZ >< GHZ| 

1/2 (|111 >< 111| + |000 >< 000| - |111 >< 000| - |000 >< 111| 



1/2 (u 3 + d 3 - 

V 2 (-^3/2,3/2,0 



~ C ) 

D (3) 
3/2,-3/2,0 



P 



(3) 



0,0,3/2 



P 



(3) 



0,0,-3/2 



For s = 0, we then have to calculate (see eq. ( 50 ) ) 

Pghz{t) 

and obtain: 



P 



GHZ 



(t) = Cl(r) P 3 ( /2, 3 /2,0 + C 2( T ) P 3/2,l/2,0 + C s( T ) P 3/2,-l/2,0 



e -3/2r e -rQ 3 g/(r)Q_ 



(3) 



P G 



HZ 



j(3) 



+c 4 (r)P, 



(3) 

3/2,-3/2,0 



-w fp (3) 

' \ / V o,o,; 



3/2 



P, 



(3) 



0,0,-3/2 



(55) 



with the weights c x (r) = 1/2 e~ 3T , c 2 (r) = 3/2 e - 2r /(r), c 3 (r) = 3/2 e - r /(r) 2 , c 4 (r) = 1/2(1 + /(r) 3 ) and c B (r) = 
l/2e _3 ^ 2T . The corresponding von Neumann entropy Sghz(t) is depicted in Figure 2. If we interpret the von 
Neumann entropy as the degree of entanglement of a state, then we see that the degree of entanglement first increases, 
and then decreases. It vanishes, once the system reaches its asymptotic state with all atoms in the ground state. 



2.5 




1 2 3 4 5 

FIG. 2: The von Neumann entropy Sghz(t) of Pghz(t) as a function of time r. 
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C. Eigenvalues of the Lindblad Operator 

In this subsection, we calculate the eigenvalues of the Lindblad operator (again with C = B/2). That is, we solve 
the equation 



LP = B[-Z/2+(l- s) fi_-(l-2a) Q 3 + s Q+]P = XP. 



(56) 



As only states with q — Z/2 have a non- vanishing trace (and as the Lindblad operator does not change q), we expand 
P in terms of states with q = Z/2: 



P 



Z/2 

£ 

q 3 = -Z/2 



r (Z) p(Z) 

L <?3 1 Z/2,q 3 ,0 



(57) 



( z) 

If we insert eq. ( 57 ) in eq. ( 56 1 , then we obtain a recursion equation for the expansion coefficients c y q3 ' . This equation 
can be brought into matrix form, which leaves us with the calculation of the eigenvalues A and the corresponding 
eigenvectors. We obtain 



X/B = 0,-1,-2,. ..,-Z 



(58) 



and the corresponding eigenvectors using the software Mathematica. If we insert the eigenvectors for X/B = 
— 1, —2, . . . , — Z into eq. (57), then we see that the corresponding states have trace zero. However, the states corre- 
sponding to the eigenvectors for X/B = have a finite trace and are therefore of direct physical interest. In this case, 
the normalized coefficients are given by the binomial distribution 



-(Z) 




,2/2+93 



(1-8) 



Z/2-q 3 



(59) 



For s = 0, the corresponding state is simply P, 



(Z) 

Z/2, -Z/2, 



= d z and for s = 1, the corresponding state is P 



z) 

Z/2, Z/2, 



For s = 1/2 the corresponding state can be compactly written as [(u + d)/2Y 



VI. CONCLUSIONS 



We have presented a group-theoretical superoperator method to solve quantum master equations for a finite number 
of atoms (or quits) . The method can be applied to many further problems in quantum optics and quantum information 
theory. Here are a few examples. For quantum optical applications, it is interesting to couple the atoms to a radiation 
field (using the Tavis-Cummings model) and to study, for example, the physics of few-atom lasers as well as phenomena 
such as resonance flourescence and optical bistability in Z-atom systems in a cavity. The resulting equations have 
to be solved numerically, but the complexity of the problems will be considerably reduced with the help of the basis 
states proposed in this paper. On the more theoretical side and with an eye on applications in quantum information 
theory, it will be interesting to construct SU (4) coherent states and study their decoherence times. Moreover, one can 
ask which many-atom states are especially stable under the influence of decoherence, and which not? We believe that 
our method will be a valuable tool in these studies. While this paper focused on fully symmetrical states, states of 
mixed symmetry and fully antisymmetrical states can also be prepared experimentally (as admixtures to generalized 
Dicke states). In a sequel to this paper, we will develop the theory of these states and study some of its applications. 



Appendix A: Proofs and commutators 



In Table II, we list the commutation relations of the superoperators Q± etc. For symmetry reasons, we list the 
commutation relations of all 18 superoperators, although only 15 of them are linearly independent. For reasons of 
space, we use the abbreviation £ = 1/2. 

To calculate the commutators [X 2 ,X' 2 ] = for all X, X' e O, wc first conclude from the commutation relations 
from Table II that 

[X 2 ,X' 2 } = [X-X + ,X'_X' + ] . 
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Q+ 


C- 




£+ 


E_ 


E 3 


M+ 


X_ 


M 3 




M- 


A/3 




«_ 




v+ 


V- 


v 3 


Q+ 





2Q3 


-Q+ 














-V+ 







U+ 







-AC 


<Q+ 





M+ 




Q- 







Q- 











V- 





iQ- 







€Q- 


N- 







-A4- 







Q 3 




















-£M_ 







-fA/l 





(M+ 







£v+ 


<V- 





s+ 











2E3 


-E+ 





U- 


-£E+ 





-V- 




-M + 







AC 







E_ 













E_ 









v+ 










M- 







-AC 


-£E_ 


e 3 

















-£M_ 







-fAC 














<V- 





M + 

















2X 3 


-M+ 














-s+ 




Q+ 





£M+ 


M- 



















A4_ 











E_ 










-Q- 


-£M- 


M 3 
































-iu. 





<v+ 


€V_ 





■A/+ 























2A/" 3 


-A/+ 


-Q+ 










s+ 


-CA/+ 


A/"- 

























N- 





S- 


-cac 


-E_ 





cac 


A/" 3 


































£v+ 


<V- 





W+ 





























2«3 













w_ 










































w 3 










































v+ 



































2V 3 


-v+ 


v_ 





































V- 


v 3 








































TABLE II: The commutators of all SU (4) superoperators. 



To calculate the remaining commutator, we first consider the case X — Q and X' — M. and obtain: 

[Q-Q+, M-M+] = M-V-Q+ - Q-V+M+ 

This expression cannot be simplified further by using commutation relations. Using the explicit definitions of the 
superoperators from eq. ([9]), we finally obtain 

[Q 2 ,M 2 ] =0. 

Similarly for the other superoperators. 
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